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Abstract 

We show how Lorentz invariance violation (LIV) can occur for Major ana neutrinos, with- 
out inducing LIV in the charged leptons via radiative corrections. Such "electrophobic" LIV 
is due to the Majorana nature of the LIV operator together with electric charge conservation. 
Being free from the strong constraints coming from the charged lepton sector, electropho- 
bic LIV can in principle be as large as current neutrino experiments permit. On the other 
hand electrophobic LIV could be naturally small if it originates from LIV in some singlet 
"right-handed neutrino" sector, and is felt in the physical left-handed neutrinos via a see- 
saw mechanism. We develop the formalism appropriate to electrophobic LIV for Majorana 
neutrinos, and discuss experimental constraints at current and future neutrino experiments. 



1 



1 Introduction 



Lorentz and CPT invariance are considered to be amongst the most sacred symmetries of elemen- 
tary particle physics. However, this very reason should motivate us to search for the smallest of 
hints of their possible violation. Indeed, motivated in part by string theory, there has been some 
recent interest in the possibility that CPT and Lorentz invariance might be violated in nature 
[1, 2]. CPT violation (CPTV) and Lorentz invariance violation (LIV) [1, 2], are clearly interesting 
effects but subject to strong constraints coming from charged fcrmions. However in the neutrino 
sector, the limits are much weaker, and so one might hope to observe such non-standard effects in 
accurate neutrino oscillation experiments [3] due to CPTV terms of the form i^fhapl'^^Li where vl 
represents left-handed neutrinos labelled by a, /? and are CPTV constants. This operator leads 
to modifications of the neutrino oscillation formula as discussed in Appendix A. A detailed discus- 
sion on Lorentz and CPT violation in the neutrino sector has recently appeared in [4], where the 
most general Lagrangian for the neutrinos in the minimal Standard Model extension is presented, 
including a catalogue of all CPTV and LIV terms [4]. 

Although of great potential interest from the point of view of future neutrino experiments, 
to stand a chance of the effects being observable, any CPTV and LIV in the neutrino sector 
must be effectively screened from the charged lepton sector, since the strong limits arising from 
charged leptons would already preclude any observation in neutrino oscillation experiments. Two 
main requirements of any effective theory of Lorentz and CPT violation in the neutrino sector are 
therefore: (i) to explain the smallness of Lorentz and CPT violation; (ii) to protect the Lorentz 
and CPT violation in the neutrino sector from the bounds coming from the charged lepton sector 
[4]. 

An elegant way of satisfying (i) is to suppose that such effects originate in the "right-handed 
neutrino" singlet sector, and are only fed down to the left-handed neutrino sector via the see-saw 
mechanism, thereby giving naturally small LIV in the left-handed neutrino sector. This possibility 
is theoretically attractive since the "right-handed neutrinos" could represent any singlet sector, 
and need not be associated with ordinary quarks and leptons, except via their Yukawa couplings 
to left-handed neutrinos. The fact that CPT violation is associated only with such a singlet sector 
could provide a natural explanation for why CPT appears to be a good symmetry for charged 
fermions, while being potentially badly broken in the neutrino sector. 

Although it is possible to satisfy (i) by appeal to the see-saw mechanism, in some cases it is 
not possible to satisfy (ii) at the same time. An example of a problematic case was discussed in [5] 
for the CPT violating operator N'^B'^^j'^N^ for the singlet right-handed neutrinos Nr labelled by 
a, (3 and the B' are CPTV constants in the "right-handed neutrino" sector. The standard see-saw 
mechanism leads to naturally suppressed CPT violation in the left-handed neutrino sector of the 
type mentioned above, namely t^lKlpI'^'^l^ where now — m\^B'^^ / [B'"^ + M^^) where rriLR is 
the Dirac neutrino mass, and Mrr is the heavy Majorana mass of the "right-handed neutrinos" . 
However the problem is [5] that the see-saw mechanism also generates unacceptably large CPT 
violation in the charged lepton sector via one-loop radiative corrections which yield the operator 
Llbioop'^f^'J^Ll, where Ll = {vl ^lY is the SU{2)l doublet that contains the left-handed charged 
leptons and neutrinos, where bgiectron ~ bf^op ^ [5]. Since the limit from the electron is 

^electron < 10~^^ eV, this implies that in the neutrino sector b < 10"^'^ eV which renders any CPT 
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Figure 1 : Tree level diagram giving rise to electrophobic see-saw suppressed LIV in the light physical left- 
handed neutrino u sector due to a lepton number violating LIV operator in the right-handed neutrino 
N sector. The dashed lines represent (possibly supersymmetric) standard model Higgs with vacuum 
expectation values < H >. The central blob represents insertion of the LIV operator in the right-handed 
neutrino sector, and the crosses indicate heavy right-handed Majorana masses. 

violation in the neutrino sector unobservable.^ 

In this paper we show how it is possible to satisfy both (i) and (ii) at the same time in a 
specific example in which the feed-down into the charged lepton sector is explicitly prevented by 
electric charge conservation. In this case it is possible to have naturally small (but still observable) 
LIV in the Majorana neutrino sector via the see-saw mechanism, without leading to any LIV in 
the charged lepton sector via radiative corrections, to all orders in perturbation theory. In order 
to ensure protection from bounds coming from the charged lepton sector the operators should 
be Majorana and lepton number violating. The essential point is that such operators, being 
lepton number violating, cannot lead to effects in the charged lepton sector due to electric charge 
conservation. As a consequence LIV could be large enough to be observable in the future neutrino 
experiments. We refer to such operators as electrophobic. We propose that such electrophobic 
operators arise exclusively from some heavy "right-handed neutrino" singlet sector, are fed down 
to physical light left-handed neutrinos via the see-saw mechanism shown in Fig. 1. The main 
purpose of this paper is to develop the formalism required for phenomenological studies of such 
electrophobic LIV operators, and to briefly study the experimental constraints at current and 
planned neutrino experiments. 

The remainder of the paper is organised as follows. In section 2 and Appendix B we argue 
that only one electrophobic operator exists and show how it can give naturally small effects in the 
left-handed neutrino sector due to the see-saw mechanism, without inducing any charged lepton 
contributions. We also derive the equation of motion for the two flavour case. In section 3 we 
derive the neutrino oscillation probabilities in the presence of electrophobic LIV, and in section 

4 we discuss the experimental constraints of electrophobic LIV at different experiments. Section 

5 concludes the paper. Appendix A contains a derivation of the equation of motion and the 

""^To be observable, the coefficient must be of the same order as an observable neutrino mass splitting y^jArr^. 
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neutrino survival probability for the usual CPT violating operator, i^L^atsl'^^L^ while Appendix B 
is dedicated to the details of the electrophobic LIV Lagrangian in a two-generation formalism. 

2 Operators, see-saw mechanism and equation of motion 

In this section we shall first write down LIV operators in some "right-handed neutrino" singlet N 
sector, arising from some high scale physics, possibly associated with the scale of heavy Majorana 
masses M. In principle M may be associated with some string scale at which LIV may be 
manifest [1], but none of our results depend on this assumption. In order to be electrophobic the 
operator should be Majoran a in nature, and violate lepton number by AL — 2. Therefore we are 
interested in operators like {N^)aOai3{Nji) f^, where are right-handed neutrinos, C represents 
charge conjugation, and O represents the remainder of the operator. There are only three such 
non-vanishing Majorana type fermion bilinears which break LIV: 



HZp{NE)aa,.{Nn)p (1) 

9%(NE)ad''a,^{Nn)3; and g'''j(M)^d-^a,,{Nn)0 (2) 

where fi, v denote the Lorcntz indices, a, ji are flavor indices and H' ^ g\ g" are LIV constants in the 
heavy "right-handed neutrino" sector. However, since the Majorana singlet that we are concerned 
with are very heavy, the terms in Eq.2 can be dropped in the static limit and the only remaining 
Majorana type LIV term which is relevant is that in Eq.l. The see-saw mechanism in Fig. 1 then 
induces LIV in the left-handed neutrino sector: ^ 

^LIV ~ Kl3{^L)a(7i,v{-I^L)l} (3) 

where h'^^ — m\j^H'^^ /{H'"^ -\- M^^^). As already noted, because of the Majorana nature of the 
operator, LIV cannot be fed down to the charged lepton sector at any loop order, due to electric 
charge conservation. By comparison the usual CPT violating operator discussed in Appendix A 
is not Majorana, and so the charged lepton sector is not protected [5]. 

In Appendix B we expand Cliv in a scenario with two neutrino flavors, and v^: 



-LIV 



(4) 



Therefore, in presence of the LIV operator that we consider, it is possible for a neutrino of a flavor 
a to transform into an antineutrino of another flavor /5, as the neutrino beam propagates. This 
would give rise to neutrino-antineutrino oscillations between different flavors due to LIV. However, 
transitions between neutrino and antineutrino of the same flavor is strictly forbidden, since CPT 
is conserved. 



^The term in Eq.3 looks similar to the magnetic moment operator for Majorana neutrinos in an electromagnetic 
field £em = fJ'a0{t^L)c^fi'^i'^L)i3F'^'' [6]. However the magnetic moment operator respects Lorentz invariance while 
the operator in Eq.3 does not. 
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The equation of motion for the two neutrino flavour case, including both the mass terms and 
the LIV terms, then follows as. 
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where 5 — l^m^ /AE, Am^ — m\ — m\ is the mass squared difference of the neutrinos and E is 
the energy of the neutrino beam. The equation of motion for Majorana neutrinos with non-zero 
transition magnetic moment in the presence of a magnetic field also has a similar form [7] . 



3 Neutrino oscillation probabilities 



In this section we look for the neutrino transition and survival probabilities in the presence of 
electrophobic LIV interactions. The neutrino mass matrix in the flavor basis can be written in 
vacuum as 

/ -(5 cos 2^ 

-5cos29 
5sin29 -H 
H 5sui29 



Mr. 
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5 sin 2^ 

-H 
5 cos 29 




(6) 



H \ 
5 sin 2^ 



5 cos 2^7 

where H is the extra element due to LIV interaction. In the above we have expressed H± — He^'^'^ 
and have assumed = 0. We will show later that is related to the change in the neutrino 
oscillation probabilities with the direction of the propagation of the neutrino and is therefore an 
important parameter [4]. However for the sake of simplicity and to get an approximate idea about 
the constraint on the electrophobic LIV term from current and planned experiments, we choose 
to put = in this and the next section. The case 7^ will be considered in section 5. The 
mass matrix in Eq. 6 can be diagonalised and the eigenvalues are 



Ai 



The corresponding mixing matrix U in the presence of the LIV term is defined as 

Mp = UUiag{Xi, A2, A3, X4)U 



(7) 



(8) 



and is given by 
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where 



Scos29 \V2 



(10) 
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One can check that when H = 0, the mixing matrix U reduces to the vacuum mixing matrix in 
the standard case and there is no mixing between the neutrino and antineutrino states. 
The general transition probabihty of a given flavor 1/^ to a flavor i/p is given by 

Pap = - 4 ^ Re{U^iU;filfip^) sin^ (^^) + 2 Im{Uo^fi;fi:fip^) sm i^\A (11) 



j>i 



where L is the distance traveled and AXij = Aj — Xj. We will assume that the mixing matrix 
U is real so that the last term in Eq.(ll) vanishes. Next we note that AA12 = AA34 = and 
AA13 = AA14 = AA23 = AA24 = 2\/6^ + lP. Thus the expression for the probability in the 
two-generation limit that we consider here reduces to 



a/3 



Sa0-4: 



UalUai + Ua2U: 



P2 



Sap - {UalUfSi + Ua2Ui, 



P2) 



sm 



(Am2)2 + {AHEf L 



We can now use Eq.(9) and (12) to get 
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(15) 



while Paa — identically, since CPT is conserved. It is again trivial to see that for if << Am^ 
or H = 0, the expressions for the probability reduces to the vacuum oscillation probabilities. On 
the other hand if if ^ Am^ then, 



P„„ = 1 - sin' ( H L 



(16) 



a/3 



(17) 



P^p = sin' ( H L 



(18) 



This implies maximal conversions of the neutrino state Ua to the antineutrino state ui^. But more 
importantly we note that the oscillations are energy independent. For the case of the usual CPT 
violating operator, ^l^plti^L^ i the survival probability given by Eq.44 in the pure CPT limit. 
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also has the same form and is energy independent [3] . Another case where the survival probability 
for the atmospheric neutrinos have the form given by Eq.l6 was considered in [9], again for a 
CPT violating theory. The form of the probability considered for LIV in [9] had a different energy 
dependence. The expressions for the survival probability that we derive here, are valid for a theory 
which does not respect Lorentz invariance, however the CPT symmetry is conserved. We derive 
the expressions for the probabilities in the massless neutrino limit, as well as for the case where 
both neutrino mass and LIV play a role in oscillations. 

The expressions Eq.13-18 are also valid for a theory with neutrino transition magnetic moment, 
in which both Lorentz invariance and CPT are conserved. However note that H corresponding to 
neutrino magnetic moment is non-zero only in the presence of an electromagnetic field. Therefore 
the case for magnetic moment is important only in the presence of an external magnetic field. 
Stringent bounds on the neutrino transition magnetic moment can be placed from solar and 
astrophysical data [7]. However the LIV term, if non-zero, is always present, irrespective of any 
other condition. 

4 Experimental constraints 

Bounds on electrophobic LIV, parametrised for example by the coefficient H discussed in the 
previous section, can be obtained from disappearance experiments using Eq.(13), and from ap- 
pearance experiments using Eqs.(14) and/or (15) ^. While the only appearance experiment with a 
positive signal is LSND, among the most prominent disappearance experiments are the solar neu- 
trino experiments, the atmospheric neutrino experiments and the reactor neutrino experiments, 
including KamLAND and CHOOZ/Palo Verde. 

Constraints from CHOOZ/Palo Verde: The CHOOZ and Palo Verde short baseline reac- 
tor experiments are consistent with no observed oscillation of Pg at baseline L ~ 1 km [10]. This 
non-observation of any oscillations can be used to constrain H^p ^ 10~^^ GcV. i/^^ (= H^p due 
to CPT invariance) is the LIV cocffccient responsible for Vei^e) '^f^i^^) transition. 

Constraints from the KamLAND experiment: KamLAND observes the electron an- 
tineutrinos produced in nuclear reactors from all over Japan and Korea. The first results from 
KamLAND show a deficit of the antineutrino flux and are consistent with oscillations [11] with 
Arn^ and mixing given by the Large Mixing Angle (LMA) solution of the solar neutrino problem 
[12]. KamLAND being a disappearance experiment is insensitive to whether the z/g oscillate into 

due to mass and mixing or z/^ due to LIV. Even though the current KamLAND data, has a 
strong evidence for suppression of the incident antineutrino flux, the evidence for energy distor- 
tion of the resultant spectrum is not very strong - no distortion of the spectrum is allowed at 
the 53% C.L. [11]. Therefore the LIV driven oscillations can explain the KamLAND data with 
H^^ ~ 7.2 X 10~^^ GeV. Though this LIV solution is not as good as oscillations with parameters 
in the LMA region, it is still allowed by the flrst results from the KamLAND experiment. It could 
be ruled out if the future KamLAND data is consistent with spectral distortion. 

^We reiterate that the bounds obtained using Eq.(13), (14) and (15) are approximate due to the neglect of the 
direction dependence of the oscillation probabilities. 
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Constraints from the atmospheric neutrino data: The atmospheric neutrino experiments 
observe a deficit of the i/^ and type neutrinos, while the observed and Pg ^-re almost consistent 
with the atmospheric flux predictions. The LIV term would convert i'fi{Pij) into Pr(^r)) while flavor 
oscillations convert v^iv^) to //^-(z/^). Since the experiments are insensitive to cither or P^, they 
will be unable to distinguish between the two cases. However, since the probability for pure LIV 
case (cf. Eq.l6) is independent of the neutrino energy, it gives the same predicted suppression 
for the sub-GeV, the multi-GeV as well as the upward muon data. This is in disagreement with 
the experimental observations. Therefore just the LIV term alone fails to explain the data and 
can only exist as a small subdominant effect along with mass driven flavor oscillations. Since the 
downward neutrinos do not show any depletion one can use Eq.l3 to put a limit of H^f ^ 10~^° 
GeV. 

Constraints from the future long baseline experiments: Better constraints on the LIV 
coefficient can be obtained in experiments which have longer baselines. The MINOS experiment 
[13] in the US and the CERN to Gran Sasso (CNGS) experiments, ICARUS and OPERA [14], 
have a baseline of about 732 km, though the energy of the z/^ beam in MINOS will be different 
from the energy of the CERN i/^ beam. However, since the LIV driven probability is independent 
of the neutrino energy, all these experiment would be expected to constrain if^^ ^ 10~^^ GeV. 
Among the next generation proposed experiments, the JPARC project in Japan [15] has a shorter 
baseline of about 300 km only, while the NuMI off-axis experiment in the US is expected to have 
a baseline not very different from that in MINOS and CNGS experiments. The best constraints 
in terrestrial experiments would come from the proposed neutrino factory experiments, using very 
high intensity neutrino beams propagating over very large distances [17]. Severe constraints, up 
to if^^ ^ 10~^^ GeV could be imposed for baselines of ~ 10, 000 km. 

Constraints from solar neutrinos: Neutrinos coming from the sun, travel over very long 
baselines ~ 1.5 x 10*^ km. So one could put stringent constraints on H^^ from the solar neutrino 
data. However the situation for solar neutrinos is complicated due to the presence of large matter 
effects in the sun. 

Constraints from supernova neutrinos: Supernova are one of the largest source of astro- 
physical neutrinos, releasing about 3 x 10^^ ergs of energy in neutrinos. The neutrinos observed 
from SN1987A, in the Large Magellanic Cloud, had traveled ~ 50 kpc to reach the earth. Neutrinos 
from a supernova in our own galactic center would travel distances ~ 10 kpc. These would produce 
large number of events in the terrestrial detectors like the Super-Kamiokande. The observed flux 
and the energy distribution of the signal can then be used to constrain the LIV coefficient. 

Constraints using the time of flight delay technique: Up to now we have been consid- 
ering the impact on the resultant neutrino signal at the detector due to spin-ffavor oscillations in 
the presence of the LIV term. The violation of Lorentz invariance could also change the speed 
of the neutrinos and hence cause delay in their time of flight. The idea is to flnd the dispersion 
relation for the neutrinos in the presence of LIV and extract their velocity v — dE/dp, where E is 
the energy and p the momentum of the neutrino beam. Then by comparing the time of flight of 
the LIV neutrinos, with particles conserving Lorentz invariance, one could in principle constrain 
the LIV coefficient. The presence of the LIV term in the Lagrangian gives a see-saw suppressed 
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correction to the mass term. Therefore 



„^ !- '»' + (19) 

where m is the usual mass of the neutrino concerned and m\jy is the LIV correction. The impact 
of the LIV correction could be important for neutrinos coming from cosmological distances. Taking 
into account the expansion of the universe, the LIV part of the mass correction introduces a time 
delay given by [18] 

At - ^ r 4^dt (20) 
2 4 a{to) ^ ' 

where a{t) is the scale factor of the universe, is the time when the neutrinos are produced, to is 
the present time and we assume that ^'q ^ 1, where 

i'l - ^ (21) 

Eq = E{to) being the energy of the neutrinos when they are observed. In principle, if one could 
estimate the At, the limit on uiliv could be used to obtain the corresponding limit on the extent 
of LIV in the neutrino sector, although, as discussed in [18], making such measurements in practice 
will be a formidable challenge. 



5 Direction dependence and the reference frame 



In this section we show that the oscillation probabilities change with the direction of propagation 
of the neutrino in the presence of the electrophobic LIV that we consider in this paper. The 
equation of motion in the flavor basis can be written in vacuum as 
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_f^g±i0_ make a co-ordinate transformation so that 



In the above we have expressed H± 
Ua — Sv'^, where 

S = diag{e"t'^'^'\ e-"^^'^/\ e"l'^'^'\ e-"^^'^/^) (23) 

Since S is a, diagonal matrix, this transformation does not change the oscillation probabilities and 
we still use the same notation for the neutrino flavor states. However the mass matrix changes to. 
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where = d(j)/dt. Thus the neutrino mixing in the presence of the LIV term we consider, and 
hence the survival and transition probabihties, will depend on 0. One can solve Eq. (24) to get 
the expression for the oscillation probabilities in presence of LIV. 

It has been stressed in [4] that in the presence of LIV interaction terms, one has to specify the 
reference frame in which the experiments are performed. They define the "Sun-centered frame" 
({61,62,63}) as standard reference frame. If we define the reference frame in which Eq. (24) is 
derived with a triad of unit vectors, p, ii and 62, where p is a unit vector along the direction 
of propagation of the neutrino and ei and 62 are the other two orthonormal vectors, then our 
reference frame is related to the standard frame through the unitary transformation [4]: 

sin © cos $ sin © sin $ cos © \ / ei \ 

cos© cos $ cos© sin $ —sin© 62 (25) 
- sin $ cos $ / V ^3 / 

where © and $ are the celestial colatitude and longitude of propagation [4]. We note that the 
angles © and $ change with the rotation of the earth and the propagation of the neutrino. This 
would make non-zero and change the oscillation probability. One can solve Eq. (24) to get the 
expressions of the mixing in presence of LIV and the oscillation probability just as we have done 
in section 3. Or one could make a co-ordinate transformation of the mass matrix given in Eq. (24) 
to the Sun-centered frame using Eq. (25) and then diagonahse it to get the oscillation probabihty 
in the Sun-centered frame. 

Thus neutrino oscillation probabilities in the presence of the electrobhobic LIV that we consider 
depend on the direction of the propagation of the neutrino. Therefore the naive bounds on the LIV 
co-efficient that we have derived in the previous section would be modified once this directional 
dependence is taken into accout. However for the most general case for i/^i, this could be quite 
an involved problem. A much more detailed discussion on the phenomenology of the i?^,^ Lorentz 
breaking terms can be found in [4]. 



6 Conclusion 

Both Lorentz and CPT violation are usually subject to very strong constraints coming from the 
charged lepton sector. Although the limits from neutrino experiments are much weaker, in some 
cases the Lorentz and CPT violation in the neutrino sector could be fed into the charged lep- 
ton sector at the one loop level, severely restricting the allowed strength of such effects in the 
neutrino sector. In this paper we have explored a class of electrophobic AL = 2 lepton number 
violating operators that induce LIV into the Majorana neutrino sector, while protecting LIV in 
the charged lepton sector to all orders of perturbation theory due to electric charge conservation. 
Among the various possible combinations, we find that the operator h^p{i'i^)a(y^y{i'L)p appears to 
be the unique candidate. This operator is Lorentz invariance violating, but it conserves CPT. To 
explain the smallncss of LIV in the neutrino sector we have assumed that LIV is introduced into 
a "right-handed neutrino" sector at some high scale, possibly close to the string scale, and feeds 
down into the left-handed sector through the see-saw mechanism, although our phenomenological 
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results are independent of this assumption. Independently of this we have developed the phe- 
nomenological formalism of the low energy electrophobic operator in the light physical neutrino 
sector. We have derived the equation of motion for neutrinos in the presence of electrophobic LIV. 
For the approximate case, where we neglect the dependence of the oscillation probabilities on 
the direction of the neutrino propagation, we have calculated the resulting neutrino survival and 
transition probabilities, and briefly discussed the constraints on electrophobic LIV arising from 
current and future experiments. We have highlighted the importance of the direction dependence 
of the oscillation probability, peculiar to the class of LIV terms considered in this paper. 
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Appendix A: The Usual CPT Violating Operator 

In this Appendix we derive the equation of motion for the previously proposed CPT violating 
operator [2]. The equation of motion in terms of the flavor states can be written as 

4("')=Hf(''') (26) 



where Hp is the Hamiltonian in the flavor basis. In this section we consider the usual CPT 
violating term considered in [2, 3], 

^2K(3^,^L^ (27) 
They argue that the only surviving CPT violating component is 6°^ (we may call it b henceforth) . 
It is a non-diagonal matrix in the flavor basis. This term has a form similar to the matter potential 
term when the neutrinos travel in matter. The Lagrangian in presence of this term has the form 

C = uH{do -a-d)u + hu^a2Mu + u^AfVau*) + u^b\ (28) 

= u\{{dQ^b^)-u ■&)u^\{u^a2Mu^u^M'^aiu*) (29) 

We note that the extra CPT violating term changes the energy component of the 4-momentum 
p^. The dispersion relation for the neutrino becomes 

£=- + 6 (30) 

where the terms have their usual meaning. We can now write Hp explicitly in this case. The 
dispersion relation (30) is actually a matrix equation and Hp is the matrix E in the flavor basis. 
We assume that and b are diagonalised by the unitary matrices Um and [/& respectively so that 

Hp = Ul^diag{ml, ml)Um + Uldiag{bi, 62)^4 (31) 
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and the Eq.(26) becomes 
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where Am^ is the mass squared difference in vacuum and A6 = 62 — ^1 is the difference between 
the eigenvalues of the matrix h. The mixing angles 9m and 6'^ correspond to the rotation angles 
that diogonalises the mass matrix and h respectively. Since there are two phases, corresponding to 
the two mixing matrices, there will be an extra phase which cannot be absorbed into the neutrino 
fields. But we have put that to zero for simplicity. It is straightforward to include it. 
We are interested in the evolution of the neutrino states. Let us define 
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Then Eq.(32) could be written as 
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That means we have two coupled differential equations 
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It is easy to solve these coupled equations using the boundary conditions 

Oee = 1 at t = 

d 

l — ttee — —A at t = 

dt 



We get 




A 



g-i(A2+B2)l/2t ^ 



A 



(^2 + ^2)1/2 

The survival probability is just the modulus squared of the amplitude a. 

A^ 



P,, = cos\A^ + B^f'h + 



(^2 + ^2)1/2^ 

3 amplitude 

sin\A^ + B^Y/H 



(^2 + ^2)1/2 

If we define 2 A = A cos 2© and 2B — A sin 2© then we get the result [3] , 



Pee = 1 - sin^ 2©sin2(AL/4) 



(33) 
(34) 



(35) 

(36) 
(37) 

(38) 
(39) 

(40) 
(41) 

(42) 

(43) 
(44) 



where A = 4((A2 + S2)V2. 
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Appendix B: The lepton number violating Cliv 

In this appendix we consider the new LIV term in Eq.3 and look for the equation of motion for 
the neutrinos. The Lagrangian for this case is, 

= C^mass + C^LIV (45) 

where Cmass contains the usual mass terms for the light neutrinos and Cliv corresponds to the 
LIV operator in Eq.3 but rewritten in 4-component Majorana notation: 



Cliv = -/i^f (*&)a(zE^-)($M)/3 + H.C. (46) 

— ^[IhtIA "^M and $m are 4-component Majorana neutrino fields. The Hermitian 
conjugate term may be absorbed into a redefinition of the coefficient as follows: 



Cliv = -<,(vI/^,).(^S^^)($m). (47) 

It is easy to see from Eq.47 that the coefficients 7Y^^ are antisymmetric and hence CPT is conserved. 
We can write this Lagrangian in the two-component notation: 

a"" = ^[a^'a" - a^a^] (48) 

^-=(T J-) (50) 

where a^^ — {I2, cr) and — — {I2, —cr). The 4-component Majorana spinor can be written in 
terms of two 2-component objects as 

*M=(|S) (51) 

where a = 1,2 and ^pa is a left-handed 2-component neutrino, while -0° = — icr^^a is the corre- 
sponding CP conjugated spinor field. For Majorana spinors 

*M = *M (52) 

And 

^ = ^mSo (53) 

= (V'd r)7o (54) 



For 



we have 



70= (J J) (55) 



*& = (r M (56) 
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Therefore the Lagrangian (47) is given by 



LIV 



n 



ab 
/if 



ab 



13^ 



(57) 



where we use m to denote the spinor with flavor a and $m to denote the spinor field of flavor h 
and have suppressed the flavor index in the 2-componcnt spinors. 

Since TYJ^^ is antisymmetric we can express it in the Lorentz space as, 



/ 


-hoi 


—ho2 


-^03 \ 


^01 





-hi2 


^13 


ho2 


hi2 





-^23 


\ho3 


-hi3 


^23 


/ 



where we have suppressed the flavor indices. Then we have 
where 



(58) 



Ho 



-Ho 



Ho — hi2 + ho3 
H± = (/i23 + hoi) ± i(/il3 + ho2) 



(59) 



(60) 
(61) 



The first term in Eq. (57) can be seen as either two incoming left-handed neutrinos of different 
flavors or alternatively as an incoming left-handed neutrino and an out-going right-handed neutrino 
of a different flavor. Therefore there is a flip of flavor as well as spin in Eq. (57). In the ultra- 
relativistic limit the full Lagrangian can be written as 



J^TJV — 



{v^r)^v,lH^ - M^i^^nH- + M^i^^rH- - (l^^^l^eLH, 



(62) 



We use this Cliv to get the equation of motion for the neutrinos. 
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